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Abstract 



1 Introduction 



The heavy ion programme at the CERN SPS was aimed at detecting signals 
for deconfined quark and gluon matter. Among the proposed signals, we shall 
consider here the enhancement of strange particle production, proposed in 
Ref. [1] and observed by NA49 and NA57, and the anomalous absorption of 
with respect to Drell-Yan pairs observed by NA50 and NA60 Collabora- 
tions, following the theoretical suggestions in Ref. [2]. 

Until now, the two signals have been separately analysed, each of them giving 
circumstancial evidence for the production of a deconfined phase of hadronic 
matter. More conclusive evidence would be provided by a correlation between 
them. Indeed, one expects that the melting of charmonia, as function of a suit- 
able critical parameter, occurs simultaneously to the reduction of the strange 
quark mass, which makes it easier for strange particles to reach thermal equi- 
librium and to appear abundantly among the hadrons formed at chemical 
freeze-out. 




Fig. 1. Strangeness undersaturation factor, 75, and J/tp anomalous suppression 
ratio, Rj/ip, as functions of the transverse size of the interaction region in heavy 
ion collisions at \/^nn = 1^.2 GeV (see text for definitions). J/ip points are from 
Pb-Pb (NA50) and In-In (NA60) collisions (triangles and circles, respectively, with 
vertical error bars only); values of 75 refer to Pb-Pb, Si-Si, C-C and p-p collisions 
(triangles, circle, box, star, with vertical and horizontal error bars), data from NA49 
and NA57. 



In this Letter we would like to point out that indeed strangeness enhancement 
and J I i/j suppression are strongly correlated when examined as function of the 
variable I introduced in Ref. [3] , related to the centrality of the collision. With 
reference to the beam axis, I is the transverse size of the overlapping region of 
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the two colliding nuclei: 
l = 2R-b, 



(1) 



R is the radius of the nuclei at rest (we consider identical nuclei) and b the 
impact parameter. 

To characterize strangeness production, we consider the strangeness under- 
saturation factor, or extra strangeness suppression factor, 75, used by many 
authors in the statistical analysis of hadron abundances. Values of 75 below 
unity mean that strange hadrons have not reached full chemical equilibrium 
in the hadron resonance gas. Thus, 75 measures how difficult is for strange 
hadrons to attain equilibrium starting from the initial collision of non-strange 
nuclear matter. 

Our main result is shown in Fig. 1, where 75 is reported as function of I, 
together with the double ratio: 



Indeed, 75 departs from its low value in proton-proton or light nuclei collisions 
to approach unity in the same centrality range where Rj/^ drops below unity. 
The possibility of such a correlation has been suggested in [5]. An analysis of 
strange particle production versus centrality was given in [6] but no correlation 
with J/ip production was considered. 

Our considerations are phenomenological in essence, however we stress that 
I is a good candidate to regulate the onset of the critical behaviour. The 
transverse lenght determines the linear size of the initial fireball and therefore 
its initial volume. It is most reasonable that collective phenomena like the 
ones we are addressing can take place only when the hadron matter involved 
reaches a certain minimum volume. In addition, the number of nucleons per 
unit surface which take part in the collision increases with /. Correspondingly, 
the energy density deposited in the fireball [4] and its initial temperature both 
increase so that, by varying /, we make a temperature scan which may lead 
into the critical or cross-over region. 



2 J/ip absorption vs. centrality 

In Fig. 2, the observed J/tp rate, normalized to Drell-Yan pairs, is reported 
as a function of I for Pb-Pb [7] and In-In [8] collisions. The J/ip attenuation 
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function is denned by [3]: 
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Fig. 2. J/ip production, normalized to Drell-Yan pairs, vs / in Pb-Pb (boxes, NA50) 
and In-In (triangles, NA60) collisions. Dashed lines represent the expected values 
after nuclear absorption correction [8]; solid lines include the hadronic absorption, 
computed assuming a limiting Hagedorn temperature: Tn ag = 180 MeV [3] (upper 
= In, lower = Pb, for both dashed and solid lines). 

p nuc i and pi are the number densities of nuclear matter and of the hadronic 
species i in the fireball, respectively; correspondingly, the cross sections <J nuc i 
and <Ti are the nucleon-J/-^ and hadron-J/?/> absorption cross sections. A ther- 
mal average is taken for the latter at temperature T [3]. The length L is the 
average linear longitudinal size of the overlapping region of the two colliding 
nuclei. The values of I as a function of the impact parameter b, for Pb and 
In, have been provided to us by the NA50 and NA60 collaborations [9] and we 
have derived the corresponding values of I with Eq.(l) and R = 6.5,5.3 fm, 
for Pb and In respectively. Finally, 3/8 • / is the average lenght traversed by a 
J/ip produced at random inside a sphere of diameter /. 

The dashed curves in Fig. 2 give the expected rates after nuclear absorption 
corrections as computed by NA60 [8]. Solid curves include the hadronic ab- 
sorption corrections, see Eq.(3), (upper = In, lower = Pb for both dashed 
and solid lines) . We outline here a few elements of the calculation of hadronic 
absorption given in Ref. [3], and refer the reader to the original paper for 
details. 

Pseudoscalar and vector mesons cross sections for: 

(ir, p, , 0, . . .) + J /if) — > open charm (4) 
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have been computed in the Constituent Quark Model [10]. The absorption 
length has been evaluated in a hadron gas at a temperature T = 175 MeV 
which was adjusted to reproduce the data points for I < 4 fm.The meson 
chemical potentials were neglected as a reasonable first-order approximation. 
Extrapolation to larger centralities requires the hadron gas Equation of State 
(EoS), to convert the energy density into temperature. We have used the EoS 
of a Hagedorn gas with limiting temperature Th = 180 MeV [3]. The ratio 
Observed (experimental points) to Expected (solid curves) is plotted vs. I in 
Fig. 1. 



3 Strangeness production 



Primary strange particle densities in phase-space are parameterized according 
to: 

p s (E,T) = (2J+l)( ls ) ns exp ( *'*~ E ) , (5) 



where n$ is the number of valence strange quarks and antiquarks in the parti- 
cle, q are a set of conserved charges, /J the corresponding chemical potentials 
and 75 is the strange particle suppression factor. In order to obtain the pa- 
rameters of the hadronic source at chemical freeze-out, i.e., temperature T, 
baryon-chemical potential jiB and 75, we have fitted full-phase-space multi- 
plicities to those measured in peripheral Pb-Pb collisions and central C-C, 
Si-Si and Pb-Pb collisions at maximum SPS energy, y/s NN = 17.2 GeV. The 
fitting procedure is described in detail elsewhere [11]; the only difference with 
respect to Ref. [11] is that we have used an updated hadronic data set [12]. A 
more detailed report on these fits is the subject of a separate publication [13]. 

For central C-C and Si-Si we have used the average multiplicities of 7r ± , K^, 
<f) and A's [6,14] measured by the experiment NA49. For central Pb-Pb the 
data sample by NA49 includes ?r ± [15], K± [15], K° [16], [17], A's [18], S's 
[19] and f2's [20]. These data samples are large enough to determine unam- 
biguously the source parameters quoted in Table 1. On the other hand, for 
peripheral collisions, the full-phase-space multiplicities measured by NA49, i.e. 
7r [21] and K [21], together with the number of wounded nucleons N w [21], 
are not sufficient to determine the source parameters. For an unambiguous 
fit, we need at least one ratio antibaryon/baryon, in order to determine the 
baryon-chemical potential. Once this is given, the ratio Nyy/n fixes the tem- 
perature and, in turn, the K/tt ratio fixes the strangeness suppression factor, 
75. Antibaryon yields as a function of centrality have been measured by the 
NA57 Collaboration [22]. These yields cover only one unit of rapidity around 
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midrapidity, however, so that some extrapolation to full-phase-space is neces- 
sary. To keep to a minimum the uncertainty related to the extrapolation, we 
have used the ratio because the rapidity distributions measured by 

NA49 in central collisions of both particle and antiparticle are gaussians to a 
very good approximation and an extrapolation is definitely more reliable than 
for A's. The procedure was as follows: 

(1) we have fitted the S + and H~ rapidity distributions measured by NA49 
in central collisions [19] to gaussians with normalization and width a left 
as free parameters and central value y = fixed; 

(2) assuming that the single gaussian shape and the ratio between their 
widths does not change as a function of centrality, we have estimated 
the yields in full phase space by rescaling the ratio at midrapid- 
ity measured by NA57 by the ratio of the widths; in fact, for a single 
gaussian, the relation between midrapidity yield dN/dy\ y=0 and N reads 
dN/dy\ y=0 = N/^hia. 

The quality of the gaussian fits is excellent and the widths of rapidity distri- 
butions turn out to be 0.860 ±0.067 for S + and 1.165 ±0.072 for S~, implying 
a ratio cr^+/a s - = 0.738 ±0.073. The obtained values of in full-phase- 

space are quoted in Table 1 for the different centrality bins. It should be 
remarked that the NA57 centrality bins corresponding to the S's yield values 
slightly different from those of NA49. However, this should imply a negligible 
error because the is fairly constant over the whole centrality range. 

The obtained source parameters for each centrality bin in Pb-Pb and for 
central C-C and Si-Si collisions are quoted in Table 1. The fit quality is 
good throughout. Temperature and baryon-chemical potential are essentially 
constant regardless of centrality and colliding system size. As expected, freeze- 
out temperatures are smaller than the temperature derived from J /if; opacity 
for / < 4 fm, which refers to the initial stage of the fireball. 

Unlike T and /i, the strangeness suppression parameter 75 varies considerably 
between the most central and the most peripheral bin in Pb-Pb. Also, 75 is 
significantly smaller in central C-C and Si-Si than in Pb-Pb collisions. These 
results are in agreement and confirm a previous analysis in Ref. [23] where a 
slightly different data set was used. 

An intriguing question about the strangeness production pattern in heavy 
ion collisions is whether a regular behaviour exists as a function of system 
size. Indeed, as can be seen in Fig. 3, 75 is larger in central light-ion colli- 
sions than in the most peripheral heavy-ion collisions, where it approaches 
the value found in pp collisions [24,13]. The search of an appropriate scaling 
variable for 75 has been recently performed e.g. in Ref. [23], where a scal- 
ing of 75 as a function of the fraction of multiply struck participants in the 
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Table 1. Particle abundances in heavy ion collisions and resulting values of thermodynamical parameters at chemical freeze-out . 
SPS energy, data from NA49 and NA57 (see text). Centrality bins are indicated for each collision. 



Glauber model was pointed out. A different approach is to look for an under- 
lying physical mechanism implying the undersaturation of the strange phase 
space, with 75 thus becoming an effective parametrization of a more compli- 
cated physical picture: for instance, in Ref. [11] a two-component model was 
proposed in which particle production stems from the superposition of a fully 
equilibrated hadron gas with 75 = 1 and single nucleon-nucleon collisions, 
where strangeness is strongly suppressed [24]; Hohne et al. [25] put forward a 
multi-cluster model where strangeness undersaturation is achieved by means 
of the so-called canonical suppression. The aim of explaining 75 in terms of a 
more fundamental physical scheme is certainly worth being pursued. However, 
it is a remarkable fact that, whatever its origin, strangeness undersaturation 
apparently scales with a pure geometrical variable such as /, as shown in Fig. 1. 




Fig. 3. The strangeness under-saturation factor, 75 as a function of the number of 
wounded nucleons in Pb-Pb (triangles), Si-Si (circle), C-C (box) and in p-p (star) 
collisions. 



4 Comparison of strangeness production and anomalous J/ip sup- 
pression 



In order to compare strangeness production to the J/ip anomalous suppres- 
sion pattern one needs to find the I values corresponding to the centrality bins 
defined by NA49. The subdivision is done on the basis of the number of pro- 
jectile spectators measured in the forward calorimeter whose relevant fractions 
[26] are quoted in Table 1. The conversion to I can be carried out by means of 
a Monte-Carlo simulation of the Glauber model. Our Monte-Carlo simulation 
of the Glauber model follows the traditional scheme [27]: at a given impact 
parameter, the number of collisions of each projectile nucleon is randomly 
extracted from a Poisson distribution whose mean is the product of the thick- 
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ness function times the inelastic nucleon-nucleon cross-section. For the target 
nucleons, the number of collisions undergone by each of them is randomly 
extracted from a multinomial distribution constrained with the total number 
of collisions as determined in the previous step and whose weights are propor- 
tional to the product of their relevant thickness function times the inelastic 
nucleon-nucleon cross-section. The thickness function has been calculated on 
the basis of a Woods-Saxon distribution: 

— = n ° (6) 

dr ~ 1 + e (r-R)/d W 



with parameters quoted in Ref. [27]: 

n = 0.17 fin" 3 ; R = 1.12A 1/3 - 0.86A^ 1/3 fm; d = 0.54 fm. 

An issue in the calculation of I is the definition of the radius R if, as it was the 
case, a Woods-Saxon distribution is used. In fact, different definitions (e.g. the 
equivalent spherical radius, square mean radius etc.) result in different values 
of I, though not dramatically. In the study of J/ip anomalous suppression, a 
value of 6.5 fm was adopted for the Pb radius, as mentioned above. To be 
consistent, we define the Carbon and Silicon radii as those encompassing the 
same fraction of the Woods-Saxon distribution as for the Pb nucleus with 6.5 
fm, namely 81% with our parametrization (6). This leads to Rc = 3.23 fm and 
i?Si = 3.84 fm. The resulting I distributions for the NA49 different centrality 
bins in Pb-Pb collisions are shown in Fig. 4. 
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Fig. 4. Probability distributions for I in the different NA49 centrality bins in Pb-Pb 
collisions, as resulting from our Glauber Monte-Carlo calculation. 
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5 Discussion and Conclusions 



Fig. 1 gives convincing evidence that neither nuclear nor hadronic absorption 
can explain the drop in J/ip production at / ~ 4 — 5 fm. In a Quark Gluon 
Plasma, the sudden drop indicated by Pb and In data would be attributed to 
the disappearance of Xc and ^'[28], which account for a sizable fraction of the 
observed J ftps in nucleon-nucleon collisions. To define precisely the expected 
signal in this case, however, one should have divided the experimental data 
by the J/vp absorption coefficient in the QGP medium. This may explain the 
residual slope displayed by the Pb points in Fig. 1, for I > 7 fm. The limiting 
hadronic absorption estimated in Ref. [3], see Eq. (3), is A -1 = E < pj(7j >~ 
0.044 fm -1 . To produce a flat behavior of the Pb points, one would need rather 
A -1 = 0.13. To our knowledge, there are no theoretical calculations available 
of the QGP absorption coefficient, which would be very useful for charmonium 
and bottomonium phenomenology at RHIC and LHC as well. 

A second remark concerns / dependence. One may wonder whether the energy 
density deposited in the fireball would be a more appropriate variable to use. 
Indeed, a plot of J/ip absorption vs. energy density estimated with the Bjorken 
formula shows better scaling of Pb and In data. Conversely, the result for 75 
looks like Fig. 3 and scaling is not achieved: the estimated energy density in 
central Si-Si collisions is in fact significantly smaller than that in peripheral 
Pb-Pb collisions, whereas the converse is true for 75 values. Errors are still 
large, however, and more precise data as well as data from other A- A collisions 
are needed for better understanding. 

In conclusion, we have pointed out that J/ip suppression and strangeness 
enhancement are correlated to the transverse size of the interaction region in 
heavy ion collisions at SPS energies, /. We cannot offer an explanation in terms 
of fundamental physics and the present result is phenomenological. However 
it is interesting to note that / determines the initial volume of the fireball 
and the energy density deposited by the collision. Collective phenomena like 
the ones we are addressing can take place most likely only when the hadron 
matter involved reaches a certain minimum volume and energy density and / 
may thus be a good candidate for a critical lenght. 

The correlation we observe makes considerably stronger the case for the SPS 
being right at the onset of QGP formation. Further experimental investiga- 
tions in this energy range are clearly called for, to elucidate the nature of the 
transition. 
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